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(D ! Abstract 

> : 

' The influence of a constant uniform external chromomagnetic fleld H on the 

5h ' formation of color superconductivity has been investigated. The consideration 

was performed in the framework of a (2+l)-dimensional Nambu-Jona-Lasinio 
model with two different four-fermionic structures responsible for < qq > and 
diquark < qq > condensates. In particular, it was shown that there exists a 
critical value He of the external chromomagnetic fleld such that at H > He a 
nonvanishing diquark condensate is dynamically created (the so-called chro- 
momagnetic catalysis effect of color superconductivity). Moreover, external 
chromomagnetic fields may in some cases enhance the diquark condensate of 
color superconductivity. 
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Typeset using REVT^ 



I. INTRODUCTION 



During the last two decades a great attention was paid to the investigation of the QCD 
ground state at finite temperature and density (see e.g. the recent review and references 
therein). The main efforts were directed to the consideration of the quark-gluon plasma 
- a new state of matter which can exist at sufficiently high temperature. Besides, it was 
also realized that at low (zero) temperature and high baryon density colored quarks, 
interacting via gluon exchange, can form Cooper pairs. Hence, the quark system would pass 
to the so-called color superconducting phase in which color symmetry of the theory is spon- 
taneously broken down. However, since the corresponding value of the diquark condensate 
< qq> was estimated to be of order 1 MeV, one could not get any observable effects in this 
case. 

Quite recently it was pointed out that due to instantons there is a nonperturbative 
mechanism of forming a condensate < qq >^ 0. As a consequence, a rather large observable 
value of order 100 MeV for the diquark condensate was predicted and the color supercon- 
ductivity (CSC) can possibly be detected in the future experiments on heavy ion collisions, 
i. e. at moderate baryon density. At present time, there is a rich literature devoted to this 
new physical effect: The CSC phenomenon has been studied in the framework of an instan- 
ton model [Q, in different versions of quark models of the Nambu-Jona-Lasinio (NJL) - 
type [P, some QCD-like theories with nonstandard color group and quark multiplets P] and 
using lattice and approaches to four-fermion models 0. CSC was also investigated in 
frameworks of renormalization group and variational as well as Dyson-Schwinger equation 
methods In all of the above cited papers [Q- the nonperturbative feature < qq >^ 
of the QCD-vacuum related to spontaneously broken chiral symmetry was taken into ac- 
count. Then, the phase structure of the theory is the consequence of a competition between 
two dynamical order parameters: < qq > and < qq >. 

It is well-known that the gluonic content of the QCD infiuences the properties of the real 
vacuum, in which there is one more nonzero nonperturbative condensate < P'^^F"-^^ > = 
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< FF >. Hence, in order to get a more adequate phase structure of the theory one should 
consider the competition of three dynamical parameters < FF >, < qq > and < qq >. 
Of course, it is very hard to solve this task within QCD itself. So, instead of this we shall 
incorporate the nonzero gluon condensate < FF > into a simpler NJL-model consideration 
of the CSC phenomenon. The NJL model does not contain dynamical gluons, hence in this 
case the gluon condensate < FF > is rather an external parameter (similar to chemical 
potential, temperature etc), than a dynamical one. In the framework of NJL models the 
condensate < FF can be realized in terms of an external (background) gauge field 

The primary goal of the present paper is just the investigation of the role which the 
gluon condensate will play in the formation of CSC. In the chosen NJL model approach we 
shall, in particular, consider a chromomagnetic gluon condensate, i.e. < FF >= > 0, 
with H being a constant chromomagnetic background field. Let us first comment the case 
of a vanishing diquark condensate. One can then imagine that the real vacuum has a color 
ferromagnet-like domain structure. Inside each domain the chromomagnetic field H"^ is 
homogeneous, but its direction is varying from one domain to another in such a way, that 
space averaging of H"^ is equal to zero. So color and Lorentz invariances are not broken [ITU 



On the other hand, inside the color superconducting phase the SUc(3) symmetry turns out to 
be broken spontaneously to SUc(2). Using pure symmetry arguments, it is easily shown that 
the three gluons living in the unbroken SUc(2) subgroup stay massless, whereas the remaining 
five gluons get masses by the Higgs mechanism. By analogy with ordinary superconductivity, 
it is expected that external chromomagnetic fields corresponding to massive gluons, i.e. 
external chromomagnetic fields of the form H"^ = {0,0,0, H^, H^), should be expelled 
from the CSC phase (Meissner effect). Moreover, sufficiently high values of such fields 
should destroy the CSC. However, our intuition tells us nothing about the action of external 
chromomagnetic fields, which in the color space look like H'^ = {H^, H^, H^, 0, 0), on the 
color superconducting state of the quark-gluon system. 

In the present paper the infiuence of such external chromomagnetic fields of the form 



H"^ = {H^, H'^, H^, 0, 0) on the phase structure of the NJL model is considered. Since we 
are mainly interested in clarifying the role of external chromomagnetic fields in the creation 
of a color diquark condensate, we put the chemical potential and temperature equal to zero. 
For simplicity, all the considerations are performed in (2+l)-dimensional space-time.Q 

The paper is organized as follows. In Section II the NJL model under consideration is 
presented and its effective potential at nonzero external chromomagnetic field is obtained 
in the one-loop approximation. This quantity contains all the information about the quark 
condensates of the theory. In the following Sections III-V the phase structure of the model 
is investigated, firstly for zero background gauge field and then for nonvanishing vector- 
potentials of two types (Abelian and non-Abelian), corresponding to the same external 
chromomagnetic field H, respectively. In Section VI it is shown that there exists a critical 
value He of the gluon condensate field at which the color diquark condensate is spontaneously 
generated (the chromomagnetic catalysis of CSC). A summary and discussion of the results 
is given in the last Section VII. Finally, detailed investigations of global minimum points of 
the effective potential are relegated to the Appendix. 

II. MODEL LAGRANGIAN AND EFFECTIVE POTENTIAL 

In the present paper the influence of a constant external chromomagnetic field on the 
phase structure of an NJL-type of model with quarks of two flavors and three colors is 
investigated at zero chemical potential and temperature. The Lagrangian of the model 
contains two different four-fermionic structures responsible for the dynamical appearance of 
< gg > as well as < gg > condensates. (Earlier, similar considerations were done for the 

^Recall that in (2+l)-dimensional space-time NJL-type models are renormalizable in the frame- 
work of nonperturbative 1/Nc expansion techniques Moreover, these theories are used in 
order to describe different planar physical phenomena, including ordinary and high temperature 
superconductivity Q|. 
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simplest NJL-type of models in which only a chiral condensate could appear 0, |T^- [0.) 
The model under consideration has the following Lagrangian: 



Here e denotes the gluon coupling constant, qc = C(f , (f = q^C are charge-conjugated 
spinors, and C = 7^ is the charge conjugation matrix (T denotes the transposition oper- 
ation). Moreover, summation over repeated indices a = 1,...,8; 6 = 1, 2, 3; /i = 0, 1, 2 is 
implied. The quark field q = qia is a flavor doublet and color triplet as well as a four- 
component Dirac spinor, where i = 1,2; a = 1,2,3. (Latin and Greek indices refer to 
flavor and color spaces, respectively; spinor indices are omitted.0) Furthermore, we use the 
notations Aa/2 for the generators of the color SUc(3) group as well as r = {t^,t^,t^) for 
Pauli matrices in the flavor space; and e and e'' are operators in the flavor and color spaces 
with matrix elements (e)*'^ = e*'^, {e^)°'^ = e"^^, where e^^ and e"'^'' are totally antisymmet- 
ric tensors. Clearly, the Lagrangian ([I|) is invariant under the color SUc(3) and the chiral 
SU(2)LxSU(2)fi groups. 

Next, let us at the moment suppose that in (|l]) the A^^^x) is an arbitrary classical gauge 
field of the color group SUc(3). (The following investigations do not require the explicit 
inclusion of the gauge field part of the Lagrangian). Below the detailed structure of A^^^x) 
corresponding to a constant chromomagnetic gluon condensate will be given. Furthermore, 
in order to demonstrate the main ideas of our approach in a way as simple as possible 
(thereby allowing for other applications in planar physics), we shall find it convenient to 
perform the following investigations in the (2+l)-dimensional space-time only.0 

Finally, notice that in order to obtain realistic estimates for masses of vector/axial- vector 



^In three dimensions the four-component spinor representation of the Lorentz group is a reducible 



one. The corresponding algebra of 7-matrices is given for example in |11,16|. Remark also that in 
our notations 7^ is an antisymmetric matrix, 7^-^=— 7^. 
^At nonzero chemical potential /x and A'J^{x) = 0, models including diquark channels similar to 
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mesons and diquarks in extended NJL-type of models [jT^, one has to allow for independent 
coupling constants Gi,G2, rather than to consider them related by a Fierz transformation 
of a current-current interaction via gluon exchange. Clearly, such a procedure does not spoil 
chiral symmetry. For the general discussion of phase transitions below (see Sections IV- VI) 
and also following the above arguments, we find it therefore convenient to treat the coupling 
constants in (|l]) as independent quantities instead of specifying them by additional model 
requirements. 

The linearized version of the model (1) with auxiliary bosonic fields has the following 
form 

L = g7^(^9^ + eA^^{x)^)q - q{a + i^'^T7^)q - + vr^) - J^A^^A'' - 

- ^*\iq^Cei'-i''q] - A''[«g£e%^Cg^]. (2) 

The Lagrangians (|l|) and (0) are equivalent, as can be seen by using the equations of motion 
for bosonic fields, from which it follows that 

a'' ~ iq^ Cee^'^^q, a ^ qq, vf ~ iq^y^fq. (3) 

Obviously, the a and vr fields are color singlets. Besides, the (bosonic) diquark field A** is a 
color antitriplet and an isospin singlet under the chiral SU(2)^xSU(2)/j group. Note further 
that a, A^ are scalars, but vr is a pseudo-scalar field. Hence, if cr 7^ 0, then chiral symmetry 
of the model is spontaneously broken. A^ 7^ indicates the dynamical breaking of both the 
electromagnetic and the color symmetries of the theory. 

In the one-fermion-loop approximation the color and chirally invariant effective action 
for the boson fields is expressed through the path integral over quark fields: 

exp(^3S'eff((T,7^,A^A*^A;:)) = J[dq][dq]exp{i J Ld^x), 



(1) were considered in low dimensions |Q] and in four-dimensional space-time as well Q, |17|. In 



contrast, we consider here the case with ^^(2;) 7^ 0, /i = 0. 



where 



5eff(cT,7^,A^A*^A;:) 



- / d^x 



a 



b /\*b 



2Gi 



+ S, 



and 



exp{i3S) = J [dq][dq] exp(i J [qDq + qMq' + q' Mq]d^x) = J [d^] exp(z J Z^d'^x). (4) 

In (^) we have used the following notations 

D = -ff'iid^ + eA^^ix)—) -a- i^^nr, M = -iA*^CeeV, M = -zA''eeVC (5) 
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and 



T /„T 



^ M , -D72^ 



V 



(6) 



L)/2 , M 

where denotes the transposed Dirac operator (see e.g. the equation (0) below). Using 
the general formula 



det 



Ik, l\ 



L , K 



det[-LL + LKL-'K] = det[- LL + LKL-'K], 



one can get from 



exp(i35) = det^/2(Z) = const ■ det^/^(L>)det^/^[D^ + AMD'^M]. 



(7) 



Let us assume that in the ground state of our model < A^ >=< A^ >=< vf >= 
and <cr>, <A'^>7^O.0 Obviously, the residual symmetry group of such a vacuum is 
SUc(2), whose generators are the first three generators of initial SUc(3). Now suppose that 
the constant external chromomagnetic field, simulating the presence of a gluon condensate 



^Clearly, < vr >^ 0, would yield spontaneous breaking of parity. In the theory of strong interac- 
tions parity is, however, a conserved quantum number, justifying our assumption that < vf >= 0. 
Nevertheless, note that at large densities a parity breaking diquark condensate could appear 
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< FF >= H'^, has the following form H" = {H\ H'^, H^,0, ...,0). Clearly, due to the 
residual SUc(2) invariance of the vacuum, one can consider the diquark condensate field 
A" = (0, 0, A3), putting = H"^ = and H'^ = H. 

Some remarks about the structure of the external condensate field y4^(x), used in (1), are 
needed. From this moment on, we select the A'J^{x) in such a form that the only nonvanish- 
ing components of the corresponding field strength tensor are ^^^2 = — -^2^ ~ i/=const. 
It is well-known that in non-Abelian gauge theories a given external chromomagnetic field 
does not fix the kind of the corresponding vector-potential uniquely. In other words, there 
can exist several physically (gauge) nonequivalent vector-potentials, which produce the same 
chromomagnetic field [T^. For example, in the case under consideration, i.e. in three dimen- 
sions, the above homogeneous chromomagnetic field can be generated by two qualitatively 
different vector-potentials: 

4(x) = (0,/ff7^,0), 4(x) = (0,0,/ff7^), A^(x) = (a = 3,..., 8), (8) 

or 

A^(x) = H6,2x'6''^ (9) 

The second of these vector-potentials defines the well-known Matinyan-Sawidy model of 
the gluon condensate in QCD pO[. In the following we shall denote expressions (§) and @ 
as vector-potentials I and II, correspondingly.^ 

There exists an attractive picture of a domain structure of the physical vacuum of QCD 
which assumes that the space is splitted into an infinite number of macroscopic domains 
each of which containing a homogeneous chromomagnetic background field generating a 



^Having background fields A^(x), given by (^) or (^), one can form three matrix fields A^{x) = 
A^^ix)^, (a = 1,2,3). Now it is easy to see that for vector-potential I (II) the corresponding fields 
A^(x) do not commute (commute) between themselves. Due to this reason, the vector-potential I 
(II) is caUed sometimes non-Abelian (Abelian) vector-potential. 
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nonzero gluon condensate < FF 0. [0. Averaging over all domains results in a zero 
background chromomagnetic field, hence color as well as Lorentz symmetries are not broken. 
(Strictly speaking, our following calculations refer to some given macroscopic domain. The 
obtained results turn out to depend on color and Lorentz invariant quantities only, and are 
independent of the concrete domain.) Note also that it was pointed out in pT| that, at high 



temperature, Abelian-like vector-potentials of the form (^ may serve as a reasonable ap- 
proximation to the true vacuum of the theory. However, at low temperature, the background 
gauge field may be essentially non-Abelian, having the form (§). 

In order to find nonvanishing condensates < a > and < >, we should calculate the 
effective potential whose global minimum point provides us with these quantities. Suppose 
that (apart from the external vector-potential A^^[x) (|)) all boson fields in S^q do not 
depend on space-time. In this case, by definition, S'efr = — Kfr / d^x, where 

V =— ln[det(D) det[D^ + AMD-^M]], v = f d^x. (10) 

Due to our assumptions on the vacuum structure, in formulae (^) we should put A^'^ = 
as well as vr = 0. Then, taking into account the relations 

Y^C = -CY, (/i = 0,l,2); X^^e^ = -e^X^, (a = 1,2,3), 

we have for the operator D (cf.(|)) with vector-potentials (|]),(0) the following identity 

D^Ce^ ^ [Y^{-id^ + eA;^(x)Aj/2) - a]Ce^ = Ce^D. (11) 

Now, using this formula as well as the relations det D = det^'y^D^'j^), det AB = det A det B, 
{ee)ij = —Sij, {€^e^)ai3 = Sa/s + ^03^/33 the expression (|T^), one can obtain after some 
evident transformations (recall that i,j = 1,2 and a,P = 1,2,3): 

V = —ln deth^ D^j^D^ + 4AA*eeeh^] 

6v 

= ^ In detl^T'^^^T'^ + 4AA*5,, (5«^ - 6^sSps)], (12) 



where A = A^. Remark, that the first term under the det-symbol in ([T2|) is a diagonal 
operator in the flavor space. One can easily see that the second term there is also a diagonal 
operator, but this time in the flavor, color, spinor as well as coordinate spaces. 



III. PHASE STRUCTURE AT ZERO EXTERNAL FIELD 

If the external field vanishes, we have the evident relation D'-j^D^^={a^ + d'^)5ij5af3- Tak- 
ing into account this formula as well as the usual one, det O = exp(tr InO), the determinant 



in (12) can be calculated straightforwardly leading to the following one- loop expression for 



the effective potential of the initial model (p!0| ) (Vcs = Vq): 



a 



2 



AA* 8i r d^k 



V„(., A. A-) = — + — + - y + 41 - k') 

where we have used the momentum space representation. This expression has ultraviolet 
divergences. Hence, we need to regularize it by introducing in (|1^) Euclidean metric — *■ 
ik^) and cutting off the range of integration (fc^ < A^). Next, by performing the integration 
and introducing renormalized coupling constants g and /, defined by 

1 4A _ 1 1 16A _ 1 

we can express the effective potential in terms of ultraviolet-finite quantities 



(14) 



Voia, A, A*) = ^ + ^ + ±(a^ + 4| AH^/^ + ^\a\\ (15) 
Zg J yTT yvr 

We shall now search for the global minimum point of the potential (|15|). Before doing 
this, let us introduce a set of more convenient parameters 

= 2|A|, A = 37T/g, fi = 37r/(2/), (16) 

in terms of which the effective potential ([T5| ) is given by 

3vrVo(a, 0) = ^ + ^ + ^(a^ + 0^)^/^ + ^1^11 (17) 
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By symmetry reasons, it is sufficient to study the function ([T7| ) in the region {cr > 0, > 0}, 
where we have the following stationarity equations 

dVo 



= a{A + 4:Ja^ + (f)^ + 2a}, 



O = 0{5 + 4J(t2 + 02}. (18) 



da 

dVo 



In order to find the global minimum point of the effective potential ([T7|) one should search 
for all the solutions of the stationarity equations (|18]) and then find among them the single 
one, where the effective potential takes an absolute minimum value. Omitting the detailed 
calculations, we present at once the results in terms of the phase portrait shown in Fig. 1. 

This figure shows the plane of parameters {A, B), divided into four domains correspond- 
ing to the four possible phases of the model. In domain 1 the point cr, = is the absolute 
minimum point of Vq, in 11 it is at values (p = 0,a = —A/6, in the region 111 the global 
minimum is located at the point cr = 0, = —B/A. Finally, in region IV the global minimum 



lies in the point a = {B - A)/2,0 = ^B^/16 - {B - A^/A. 

Recall, that the coordinates of the global minimum point of the effective potential are 
the vacuum expectation values < o" >,<</>> of the fields a, (p. Lagrangian @ provides us 
with the equations of motion for a, A from which one can easily get the following relations 
< cr >~< qq >, < (f) >~< qq >, i.e. the global minimum point of the effective potential 
gives us an information about chiral and Cooper-type diquark condensates of the model. 
Hence, if the parameters {A, B) belong to the region 1, we have the symmetric phase, because 
of < cr >=< qq >= 0, < (j) >=< qq >= in this case. The phase with spontaneously broken 
chiral symmetry is situated in region 11, since here the chiral condensate < qq >^ 0. In 
region 111 the diquark condensate < gg > is nonzero, thus indicating the presence of the color 
superconductivity phase. Clearly, in this case the electromagnetic as well as color symmetries 
of the model are spontaneously broken. Finally, region IV corresponds to the mixed phase 
of the theory, where chiral, electromagnetic and color symmetries are spontaneously broken 
down (here < qq > , < qq >^ 0). 
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IV. PHASE STRUCTURE OF THE MODEL FOR VECTOR-POTENTIAL I 



Let us next study the influence of the external chromomagnetic field with constant vector- 
potential @ on the phase structure of the model (1). By using the momentum space 
representation for the operator under the det-symbol in (|12|), we obtain instead of the 
differential operator an algebraic (24x24) matrix which has three different eigenvalues Ei(p) 
(i=l,2,3) 

eH , 1 



Ei,2{p) = r - PS + ^ + 4| + ^ ± ^\/ieHy + AeHp^, E,{p) =p'-pi + a\ (19) 

each of them having an eight-fold degeneracy, and = p\+p\- Taking into account this fact, 
one can easily obtain from ( [10|) and ([T2|) the following expression for the effective potential 
of the model in the presence of an external vector-potential (H) (Kg = Vh^)\ 

^2 AA* AiJ^ f (Pp 



a 



In E,{p), 



(20) 



2Gi ■ G2 ■ 3 (2vr)3 

Integrating here first over pq and then over the space of variables and employing a 
suitable ultraviolet cutoff, one can get after adopting a renormalization procedure (similar 



calculations were performed in [|T^ for the model (1) in the case with G2 = 0)'- 



SttVh, {a 
3eH 



Aa^ B4)^ 2 a P 2 r 

\ — + + - 

2 2 3 3 



Cr 



(T^ + 0^ + eH) 



1/2 



(21) 



where we have used the notations introduced in formula (P!B|). Instead of the dimensional 
quantity (PT|), let us consider the dimensionless function Vi{x,y) = 37rVf/-^(cr, 0)/(eif)'^/^, 
where x = a/ VeH, y = (j)/ V eH. Evidently, we have 



Vi{x,y) 



Ax^ 



+ 



By 



2 r 



2 3' ' 3 



^;i + x2 + y2__(a;^ + ^^)ln 



1 + Vl + x^ + y^' 



(22) 



where A = A/V^, B = B/V^. 
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Phase structure in terms of A and B 

In contrast to the potential (0), which has three dimensional parameters A, B, eH, the 
function ( p2D depends only on two parameters A, B. So, first of all we shall study the 
global minimum dependence of the potential Vi on the parameters A, B. Doing this, we get 
the possibility to discuss the phase structure of the model in the presence of the external 
vector-potential (|]). Since Vi is symmetric under transformations x —x or y ^ —y, it 
is sufficient to look for its global minimum only in the region x,y > 0. The stationarity 
equations then take the form 

—^ = x{A + 2x + Fiz) I r^-^ } = , — i = y{ B + Fiz) I r^-^ } = 0, (23) 

where 

F{z) = 2z + 2Vl + ^2 _ (1 + Vl + z'^)/z ]. (24) 

We should find all the solutions of the equations (^3]) in the region x, ?/ > and then select 
that one, where the function Vi{x,y) takes its global minimum. Omitting here calcula- 
tional details, we directly quote the results in the form of the phase portrait. The detailed 
investigation is given in Appendix A. 

In Fig. 2 the phase portrait of the model (1) in the presence of a nonzero vector-potential 
of the type I is presented in terms of A and B. Here one can see three regions. Above the 
line li there is the Phase III of the model, which corresponds to the Vi (x, y) global minimum 
point of the form (0,yo(-S))- (The properties of the function yo{B) as well as the functions 
xo{A), xi{A, B), yi{A, B) considered below are given in Appendix A.) In this case < a >= 0, 
< (p >7^ 0. Below the curve I2 the global minimum of the effective potential has the form 
{xq{A),0). Thus, in this region the Phase II is located, since for such values of A,B the 
model has a vacuum with < a >^ 0, < >= 0. Finally, inside the i7-domain there is 
a mixed Phase IV, since here the global minimum point {xi{A, B),yi{A, B)) corresponds 
to the vacuum with < qq and < qq 0. It is necessary to emphasize that in the 
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presence of such kind of external vector-potentials the Phase I is absent at all (obviously, 
only when ^^(a;) = 0, this phase can be realized in the model). 



V. PHASE STRUCTURE OF THE MODEL FOR VECTOR-POTENTIAL II 

Next, let us study the influence of a nonvanishing external chromomagnetic field with 
vector-potential on the phase structure of the model (1). In this case, after some calcu- 
lations, the operator, which stands under the det-symbol in (|12|), can be transformed to the 
following expression in the color space: 

D7^D7^ + 4AA*5,,(5„^ - 6^,6p,) = diag(5,,D+ , 6,,D^ , 6,,{a^ + d')), (25) 

where 

IP— — — 

D± = a^ + 4| A|2 - n±n±^ T jYYF,., nj = td, ± eA^{x)/2, = HS^^Xi- (26) 



Remark, that in (|2|), (|6D D± are operators in the coordinate and spinor spaces, only. The 
same is true for the expression (a^ + d"^), which is the unit operator in the spinor space. 



Taking into account (^), one can easily find for the potential V (|T2D the expression: 



V = — tr In D+ + — tr In + — tr ln(cr^ + d^) (27) 

3v 3v 3v 

(the trace prescription in (p7|) is taken over coordinate as well as spinor spaces). The 
last term in this formula was calculated in Section III of the present paper. Concerning 
the first two terms in (^), we should point out, that D± = D±'~f^D±'~f^, where D± = 
ij'^d^ ± e'~f^A^{x)/2 + M are formally the Dirac operators for fermi-particles with electric 



charges ±e/2 and effective mass M = ya'^ + A\A\'^ in the presence of a constant external 
magnetic field H. Similar expressions were calculated in a lot of papers (see, e.g. |2^, from 



which it follows that the first term in (|27|) equals to the second one). So, we omit details of 

tr In D± calculations and present at once the corresponding effective potential of the model: 

(T^ AA* pH r°° fl^ 
V'„,KA,A-) = — + — + ^/^ _exp(-.(a' + 4|An) coth(eHs/2) 
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In this formula eH has a positive value. Both integrals in (^8]) are ultra-violet divergent 



ones. To renormalize the first integral one can act in the same manner, as it was done in |T6 
with the effective potential of the three-dimensional Gross-Neveu model in the presence of 
an external magnetic field. The second integral in ( pS]) was already renormalized in Section 
III. Hence, the finite expression for the effective potential of the model (1) in an external 
chromomagnetic field of type II looks like: 



VH.(cr,A,A*) = Vo(fx,A,A*) + 
eH r°° ds 



coth 



feHs\ 2 



V 2 V eHs. 

where Vo(cr, A,A*) denotes the effective potential at if = (see (|15D). Integrating in 



(29) 



this formula with the help of integral tables [^], one can get the following more compact 
expression for the effective potential 

(T^ AA* 2 



VhM a, a*) = - + ^ + -i^r + 

2g f 9n 



eH^a^ + A\A\^ 2{eHfl\ ( 1 a^ + 4|Ap \ 
+ Stt Stt 4 2' eii ' ^^"^ 



where C,{v^x) is the generalized Riemann zeta-function ^4[. As in the previous section, let 
us further introduce the dimensionless function V2{x,y)= 3TT(eH)^^/'^VH2{o', A, /S.*): 

V2{x, y) = ^ + ^ + '^\x\' + v/x2 + y2 _ 2C(-l/2, + y^), (31) 

where x = cr/V eH, y = 2| A|/V eH and A,B are the same parameters, as in the relation 



(Ua)- Clearly, this expression differs from the corresponding quantity (|2^) for a non-abelian 
background field. 

The reader, who is not interested in following the details of our investigation of a global 
minimum point for V2{x, y), can at once look at the phase portrait of the model in terms of 
A, B. Qualitatively it is the same as the phase portrait for the function Vi{x, y) (see Fig. 2), 
i.e. it contains only three different Phases II, III and IV. Details of calculations are quoted 
in Appendix B. 
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VI. CHROMOMAGNETIC CATALYSIS OF COLOR SUPERCONDUCTIVITY 



Let us now analyse in more details the phase portrait of the model (1), this time in terms 
of A, B, eH. In particular, we shall describe phase transitions which occur for arbitrary 
fixed A, B and with varying values of H . In general, our discussions concern both cases 
with vector-potentials of types I and II simultaneously. However, where it is necessary, we 
indicate to what type of vector-potential the information is related. First of all a general 
remark: If A, B are fixed and H is varied from to cxd, then in the plane (A, B) one moves 
along some ray (which depends on A, B) from infinity to the origin (this fact simply follows 
from the definition of A, B). 

The case A, i? > 0. (In this case we have a weak coupling for both bare constants: 
^1,2 < Cc ~ TT^/A.) At = this choice of parameters corresponds to the unbroken 
Phase I (see Fig.l). If A, B are fixed in such a way, that A> B {in terms of bare coupling 
constants this means ^ > + ^)-, then at if 7^ we have in the {A^B) plane of Fig. 
2 a ray which is located above the line /i, i.e. is in the Phase III. Hence, in this case the 
external chromomagnetic field induces (catalyses) the dynamical generation of a nonzero 
diquark condensate. Here in the point He = 0+ one has a second order phase transition 
from Phase I to the phase with color superconductivity. At varying values of H the diquark 
condensate behaves, e.g. in the case with vector-potential II, in the following wayQ: 

< qq >~ feH at H ^0 , <qq>^ VeH at H ^ 00. (32) 

(The chiral condensate in this case is identically zero.) 



^In order to find < qq >^< A > and < qq >r^< a >, one should multiply the coor- 
dinates of the global minimum point of the functions Vi{x,y) or V2{x,y) (see Appendixes A 
and B, correspondingly) by the quantity \/eH. For example, in the case under consideration 
< qq >~< A >= yQ{B)\/ eH/2, where the function yo{B) and some of its properties are presented 
in formula (B.4) of Appendix B. 
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U A < B and the external chromomagnetic field H varies in the interval H e (0, oo), 
then points in the {A, B) plane vary along a ray r (see Fig. 2). If if ^ 0+, we are at the 
infinite end of this ray, i.e. in the Phase II of the model (see Fig. 2). Hence, in the point 
He — 0+ the chromomagnetic field induces a dynamical chiral symmetry breaking phase 
transition (a second order phase transition), but a diquark condensate is not produced. 
These properties of the vacuum are not changed for sufficiently small values of H such that 
H < Hc{A, B). The value H — Hc{A, B) corresponds to the point o (see Fig. 2) in which the 
ray r crosses the line I2 and passes to the Q-region where chiral and diquark condensates are 
both nonvanishing. So, with growing values of H in some point Hc{A, B) one has a second 
order phase transition from the Phase II to the mixed Phase IV. In Table 1 the results of a 
numerical investigation of H(.{A, B) as a function of A are presented for some fixed values 
of B/A in the case of a vector-potential of type I. One can see from this table that if the 
ratio B/A is fixed, then Hc{A, B) increases with A as A^. It is also clear that for each fixed 
value of A the quantity Hc{A, B) is a growing function of B. 

The behaviour of condensates in the case A < B and for a vector-potential II are the 
following: 

< qq >~ geH at ^ 0, < >~< a >= ~ at H^{A, B) < H, (33) 

<qq>=Q at H < Hc{A, B), < qq >^ \/eH at H ^ 00. (34) 

So, at H ^ the Phase I of the model is completely absent in the phase structure of the 
model for both types of external chromomagnetic fields. 

The case .4 < 0. 2.4 < 3D. In this case at H — one has a Phase II of the theory (see 
Fig. 1) with spontaneously broken chiral symmetry. If the external chromomagnetic field 
H varies in the interval H e (0, 00), then in the {A, B) plane there is a ray which crosses 
the line I2 in some definite point. If if — > 0+, we are in the infinite end of this ray, i.e. in 
the Phase II of the model. However, in contrast to the previous case with A,B>Q and 
A < B/va. the present case the value if = 0+ is no more the point of a phase transition. (At 
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A < 0, when the bare couphng constant Gi takes a supercritical value Gi > Gc, the origin of 
chiral symmetry breaking is the rather strong supercritical quark- ant iquark attraction, but 
not the chromomagnetic field. In this case the external chromomagnetic field only stabilizes 



the vacuum with chiral symmetry breaking [0,|TB].) If H increases, we move along this ray 



to the origin of the (A, i?)-plane. Hence, starting from some value Hc{A,B), we are in the 

region Q (see Fig. 2), where besides < qq >^ the diquark condensate is nonzero as well. 

So, at sufficiently high values of H > Hc{A, B) the Phase II of the theory is transformed 

into a mixed Phase IV. 

The infiuence of the vector-potential II on the chiral Phase II of the model is realized in 

the following behaviour of condensates: 

A (B — A] 

< qq >~< o- >= -— at IT ^ 0, < qq >~< a >= - — - — - at Hc{A, B) < H, (35) 
6 2 

<qq>=0 at H < Hc{A, B), < qq >~ VeH at H ^ oo. (36) 

The case B < 0, A > B. In this case at H = there is a perfect (not mixed) color 
superconducting Phase III of the theory (see Fig. 1). One can easily show that now for 
all values of H only the diquark condensate < gg > is nonzero. This vacuum is chirally 
invariant, but the Uem(l) as well as color SUc(3) symmetries are broken down. It is possible 
to show that in this case 

< qq >~< A >= -B/4: at H ^0 , < qq >^ \feE, at H ^ oo, (37) 

i.e. the external chromomagnetic field even enhances the color superconductivity. 

The case A < B, 2A > 3B. Analyzing here the behavior of the quark condensates in a 
similar way as in the previous cases, one can easily establish that the vacuum properties 
are not changed with growing values of H. Hence, at if = as well as at if 7^ there is 
a mixed Phase IV with nonzero quark and diquark condensates. Remark that the action 
of an external chromomagnetic field on the mixed phase does not change the value of the 
chiral condensate: It is the same as at ii = 0, where < qq >~< a >= {B — A) /2. However, 
the diquark condensate depends on the value of H: 



< qq >~< A >= Y^^Vie -{B- A)2/4 at if ^ 0, < qq >~ VeH, at H ^ oo (38) 

In conclusion, let us remark that for arbitrary fixed parameters A, B and in the presence 
of sufficiently large values of external chromomagnetic fields of both types there arises a 
nonzero diquark condensate < qq >^ 0, i.e. the color superconducting phase of the model is 
realized. \i A> 3i?/2, i? < (i.e. for sufficiently high values of G2 > Gc), then < qq >^ 
even at ii = (in this case the external chromomagnetic field enhances the CSC). However, 
for other regions of the (A, i?)-plane the nonzero external chromomagnetic fields catalyse 
the generation of < qq >^ 0. The critical value of if, at which color superconductivity is 
induced, may be 0+ (if A, i? > 0, A> B), or some finite value Hc{A,B) 7^ (in the last 
case we have not a perfect, but mixed color superconducting phase in which the diquark 
condensate coexists with the chiral condensate < qq >^ 0). Apart from this, in the presence 
of a chromomagnetic field, the phase portrait of the model does not contain the symmetric 
phase. 

VII. SUMMARY AND CONCLUSIONS 

In the present paper the phase structure of a (2+l)-dimensional four-fermionic NJL-type 
of model (1) with two coupling constants was investigated admitting nonzero background 
vector-potentials of two nonequivalent types I and II (see (|])-(^). In the framework of such 
a model the external vector-potential might be thought to simulate such a nonperturbative 
feature of the real Q CD- vacuum like a nonzero gluon condensate < FF >= H^. The 
structure of the Lagrangian (1) permits us, in particular, to consider then the competition 
between chiral < qq > and diquark < qq > condensates and to get some insight into the 
role of the gluon condensate as a possible catalyst of color superconductivity. 

It is well-known that color-superconducting quark matter with two quark flavors arises 
by the condensation of color antitriplet diquark Cooper pairs. The condensate breaks the 
SUc(3) symmetry down to SUc(2). Hence, the three gluon fields corresponding to the gen- 
erators of unbroken SUc(2) stay massless and the remaining five gluon fields receive a mass 
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by the Higgs mechanism (Meissner effect). We have studied the inffuence of the exter- 
nal chromomagnetic fields living in an unbroken SUc(2) subgroup of SUc(3), i.e. having a 
form = {H^ , H^, H^,0, ...,0), on the formation of the color diquark condensate. Using 
a global SUc(2) color rotations one can bring this field to the form if" = (0, 0, if, 0, 0) 
which corresponds to the above mentioned vector-potentials (§|)-(^). 

The main conclusion from our investigations is that at zero chemical potential the ex- 
ternal chromomagnetic fields of these type are a good catalyst of color superconductiv- 
ity. (Earlier, it was shown that external (chromo) magnetic fields catalyse dynamically 
the spontaneous breaking of chiral symmetry in some (2+l)-dimensional four-fermionic 
models p2| , p5| , p!3| , p!6| . It turns out that this is a particular manifestation of the so-called 
magnetic catalysis effect (see e.g. [p6| , p!5| , p7| ), which has a rather universal model indepen- 



dent character.) Indeed, we have shown that for sufficiently small bare coupling constants 
Gi^2 < TT^/A, i.e. for such values of Gi^2 at which for H = one has a symmetric Phase 
I of the theory (see Fig. 1), the pure CSC phase (< qq >= 0, < qq >^ 0) is realized in 
the model at infinitesimally small values of the external chromomagnetic field H ii G2 > Gi 
(in terms of A,B that means A > B). If G2 < Gi {A < B), then, first, a chirally breaking 
phase transition induced at if = 0+ (chromomagnetic catalysis of chiral symmetry breaking 
at which < qq >^ 0, < qq >= 0) occurs. After that, with growing values of H, at some 
point H = He there is a second phase transition to the mixed phase of the theory, where 
both condensates < qq > and < qq > are nonzero (both phase transitions are continuous 
second order ones). 

The action of an external chromomagnetic field on the chiral Phase II of the theory (see 
Fig. 1) is to induce < qq >^ at some critical point H = He 0, thus drastically changing 
the vacuum properties and transposing the system into a mixed Phase IV. 

Finally, we should mention that the ground states of Phases III and IV are not changed 
under the influence of an above mentioned external chromomagnetic field. So, the external 
chromomagnetic fields living in an unbroken SUc(2) subgroup of SUc(3) only enhances the 
CSC phenomenon. 
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Notice that all the above mentioned effects are observed in the presence of both vector- 
potentials I, II. 

To our opinion there exists a deep connection between the chromomagnetic catalysis of 
color superconductivity and chiral symmetry breaking, induced by external chromomagnetic 
fields ||TB[- ||TB[. This assumption is based on the existence of the Pauli-Giirsey (PG) transfor- 



mation ||2^ , mixing quarks and antiquarks, due to which some phenomena in the gg-channel 
can have its analogy in the gg-channel. In particularly, this suggests that diquark conden- 
sation might be understood as the properly PG-transformed chiral condensation. However, 
the detailed consideration of this question is not the subject of the present paper and will 
be investigated elsewhere. 

Moreover, in the nearest future we are going to include into our consideration of the 
simple NJL model (1) a nonzero chemical potential /i in addition to the external chromo- 
magnetic fields. Recently, in the framework of NJL models ||2^ the influence of /i and an 
external magnetic field on the chiral properties of the vacuum were considered. Apart from 
discovering different kinds of magnetic oscillations (relativistic van Alphen-de Haas effect) 
in the strongly interacting quark systems, there it was also found that in the NJL model at 
nonzero baryon density the chiral symmetry must be restored at sufficiently large values of 
the magnetic field. By analogy with we expect that the diquark condensate < qq > 
should disappear in the case of nonzero baryon density for a sufficiently strong external 
chromomagnetic field, i.e. at large values of the gluon condensate. 



ACKNOWLEGMENTS 

We wish to thank V. A. Miransky and Y. Nambu for interesting discussions and A.K. 
Klimenko for numerical calculations. One of us (D. E.) acknowledges the support provided 
to him by the Ministry of Education of Japan (Monbusho) for his work at RCNP of Osaka 
University. This work is supported in part by the Russian Foundation of Basic Researches 
(project 98-02-16690) as well as by DFG-Project 436 RUS 113/477/4. 



21 



Appendix A 



Investigation of the global minimum point of Vi{x,y) 

It follows from in the text that the function F{z) monotonically increases on the 
interval z G (0, oo) and F{0^) = — oo, F{+oo) = +00. Hence, for arbitrary fixed values 
of A,B there exist only two real numbers Xo{A) > 0, yo{B) > 0, such that the two pairs 
(xo,0), (0,?/o) as well as the trivial one (0,0) are solutions for the system of stationarity 
equations p3|). (The Xo{A) and yo{B) are zeros of the functions which are located inside the 
first (at y = 0) and second (at x = 0) pair of braces in (p3[), respectively.) Furthermore, one 
can easily see that dVi/dx < if y = and x G (0, Xo{A)) as well as dVi/dy < if x = 
and y G (0,?/o(-B)). This means that the quantities Vi(a;o,0) and Vi(0,?/o) are smaller, than 
Vi(0,0). So, for arbitrary finite real values A,B the global minimum point of the function 
Vi{x,y) cannot lie in (0,0). Due to this reason the symmetric phase is absent at all in the 
phase structure of the model. 

In the next formulae some properties of Xq{A) and yo{B) are presented: 

xo{A) = 2e-^-2 at i ^ +00, so(0) = 0.147..., xo{A) = -A/6 at i ^ -00, [A.l) 

yo{B)=2e-^-^ at B ^ +00, yo(0) = 0.183..., yo{B) = -B/4: at B ^ -00. (A2) 

It follows from (|23|) that there may exist (but not for all values of A, B) one more solution 
{xi{A, B) ,yi{A, B)) of the stationarity equations, where xi > 0,yi > 0. (For this solution 
the functions which are inside both braces in (^) take zero values.) Evidently, we have 

xM, B) = {B- A) II, |/i(i, B) = ^yl{B) - {B - Af/A. (A3) 

From (A. 3) one can easily see that this type of solutions for equations ( P5D exists inside the 
region Vt of the (A, i?)-plane (see also Fig. 2): 

Vt = {{A,B) ■ B -2yo{B) < A< B}. (AA) 
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Using equations ( ]23|) one can find the following values of the potential (p2D in its stationary 
points: 

V,{x,{Al 0) = _ 2xl{A)) - \xl{Al (A.5) 

Vi{0,yom = )ll±^^ii _ 2yl{B)) - \yl{Bl (A.6) 

V,{xM^B)MA.B)) = _ 2yliB)) - ^yUB) - (A7) 

On the line h = {{A, B) : A = B}, which is a part of the boundary for the region Q, we have 
= (0,?/o)- Hence, on this line Vi(0,yo) = Vii^ijUi)- Comparing (A. 6) and (A. 7), we 
see that inside the f2-region Vi(0,?/o) >^i(a^ijZ/i)- 

The other part of boundary for the region Q is the line I2 = {{A, B) : A = B — 2yo(-B)}. 
With the help of the stationarity equations it is possible to show that on this line the 
following relations are also fulfilled: B = A + 2xo{A), xo{A) = yo{B). As a consequence, we 
have {xi,yi) = (xcO) as well as Vi(a;o,0) = Vi(xi,yi) on the line l2- 

Numerical investigations show that inside the region Q there is a line on which 
Vi(xo,0)= Vi{0,yo). Further, it is important to remark that the derivative of the func- 
tion Vi(0, ?/o(-S)) with respect to B as well as the corresponding partial derivative of the 
function Vi{xi{A, B),yi{A, B)) are positively defined quantities in their regions of defini- 
tion. Now, taking into account all the above mentioned facts, it is possible to assert that in 
Fig. 2 the phase portrait of the model (1) in the presence of a nonzero vector-potential 
is presented in terms of A and B. This means that above the line /i there is the Phase III of 
the model, which corresponds to the Vi{x,y) global minimum point of the form (0,?/o)- (Iii 
this case < a >= 0, < >^ 0.) Below the curve I2 the effective potential global minimum 
has the form (xq, 0). So in this region the Phase II is located, since for such values of A, B 
the model has a vacuum with < a >^ 0, < >= 0. Finally, inside the f2-domain there is 
a mixed Phase IV, since here the global minimum point {xi,yi) corresponds to the vacuum 
with both nonzero condensates < qq > and < qq >. 
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Appendix B 



Investigation of the global minimum point of V2{x,y) 

In this Appendix we present the search of the global minimum point for the potential 
V2(x, y) as well as its dependence on the parameters A, B. Since this function is symmetric 
under two discrete transformations x — *■ —x and y — > — y, it is sufficient to study it only in 
the region x, y > 0. The stationarity equations for ^(a;,?/) take the form 

dV2{x, y)/dx = x{A + 2x + (x^ + y'^)"^''^ - 2C(l/2, x^ + y^)} = 0, (5.1) 

y)/dy ^ y{B + (x^ + y^Y^I^ - 2C(l/2, x^ + y^)} = 0. {B.2) 

One can see from (B.1)-(B.2) that the first derivatives of V2 do not exist in the point (0, 0). 
(This means that if the point (x, y) tends to the origin along different ways, the resulting 
expressions for the partial derivatives at the point (0, 0) do not coincide.) In contrast, the 
function Vi{x,y) is differentiable in the point (0,0). So, we need a special investigation of 
this point. Let us put ?/ = in the equation (B.l). Then, using properties of the (^{I'.x)- 
function it is easily seen that at ?/ = and x ^ 0+ the partial derivative dV^jdx tends 



to (-1). Analogously, at x = and ?/ ^ 0+ the derivative dV^jdy (B.2) tends to (-1) as 
well. This fact means that for arbitrary values of A, B the point (0, 0) cannot be a global 
minimum for the potential V2(x,?/). So, in contrast to the case with H = 0, the ground 
state with intact initial symmetry is no more possible in the model (1) at H ^ 0. Such 
a property of the effective potential is a characteristic feature for a phenomenon which is 
called (chromo)magnetic catalysis of dynamical symmetry breaking. According to this effect 
the external (chromo)magnetic field promotes in a great extent the spontaneous breaking of 
initial sjTiimetry of the theory (for more details, see the last section of the present paper). 

Similar to the case with non-abelian vector-potential of the type I, in the present con- 
sideration it is possible to show that for arbitrary values of A, B the stationarity equations 
(B.1)-(B.2) have two solutions of the form (xo(A),0) and {0,yo{B)), where 

xo{A) = l/A at xo(i) = -i/6 at i ^ -00, (5.3) 
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yo{B) ^l/B at B ^ +00, yo(^) = -S/4 at S ^ -00. (BA) 

(Here and in the following discussions of the present Appendix we use the same notations 
xq{A) and yo{B) for the solutions of stationarity equations as in the previous Appendix. 
But one should remember that these functions have quite different numerical values, than 
similar functions had in Appendix A.) 

Prom (B.1)-(B.2) it follows that only for {A, B) e where Q is defined formally in (A.4), 
there is a solution of the form {xi{A, B),yi{A, B)), where Xi > 0,yi > 0. These functions 
are given in (A. 3). There are no other solutions of the stationarity equations. 

Using numerical and analytical methods, it is now possible to compare the values of the 
effective potential in its stationary points and thus to find the global minimum point of V2 
as well as its dependence on the parameters A, B of the theory. We omit the details of 
this investigation and present only the results, which can be formulated in form of a phase 
portrait of the model in terms of A, B. 

It turns out, that Figure 2, which is the phase portrait of the model for nonzero vector- 
potential I, formally may serve as a phase portrait of the model for nonzero external gauge 
field of the type II as well. In both cases the line h separates the color superconducting 
Phase III from the mixed Phase IV. Further, in both cases, the Phase IV is separated from 
the chiral Phase II by the line I2, which has the same analytic definition through xo{A) and 
yo{B) (see the figure caption to Fig. 2). Moreover, the leading asymptotic behaviours of the 
/2-curves are in both cases identical (at yl — > 00 we have instead the line /i, at A ^ —00 it 
is the fine 2 A — SB). Of course, since for type I and type II vector-potentials the functions 
Xo{A) and yo{B) obey different stationarity equations, the line I2 of case I does not coincide 
with the line I2 of case II. Finally, we should stress again that for both types of nonzero 
vector-potentials I and II, the symmetric Phase I of the theory, which is present in the phase 
structure of the model at if = (see Fig. 1), is absent at all. 
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Figure captions 



Fig.l Phase portrait of the model at vanishing chromomagnetic field H — 0. The 
boundaries h, h are defined by h = {{A, B): B}, h = {{A, B) : 2A^ W}. 

Fig. 2 Phase portrait of the model in the presence of a nonzero vector-potential 1. The 
symmetric Phase I is absent at all. The boundaries li and I2 of the region Q are defined 
according to [i = {{A,B) : A ^ B}, h = {{A,B) : xo{A) = yo{B)}. At B ^ +00 
[B —00) the curve I2 approaches asymptotically the line /i (the line 2 A = 3B). The 
curve I2 intersects the axes A and B in the points (-0.36) and (0.3), respectively. For fixed 
A and B and varying H one moves in the {A, S)-plane along some ray r. At if = Hc{A, B) 
there is a phase transition from the Phase 11 to the mixed Phase IV, if the ray r intersects 
the line I2 in some point o (see the detailed discussion in the Section 6). 



29 



TABLES 

TABLE I. The ^-dependence of Hc{A, B) for some fixed values of the ratio B/A. 



B/A 


100 


10 


2 


1.5 


1.1 


1.05 


eH,iA,B)/A^ 


112925.63 


967.55 


15.804 


5.1038 


0.6593 


0.3508 
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